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numerical results for the exact amplitude. 
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1. Introduction 

The top quark, with a mass of approximately 173 GeV, is the heaviest elementary particle 
produced at colliders until now. So far, the study of the properties of the top-quark was only 
possible at the Tevatron, where the mass of this particle was measured with an accuracy 
of less than one percent. The production cross sections and decay widths are only known 
with larger uncertainties. Because of its large mass, the top quark is expected to couple 
strongly with the electroweak symmetry breaking sector. Therefore, the study of scattering 
and decay processes involving top quarks is expected to provide fundamental clues on the 
mechanism responsible for the origin of particle masses. This is one of the main goals of 
the Large Hadron Collider (LHC) physics program. 

While the Tevatron produced a few thousand top quarks, the LHC will be an authentic 
top-quark factory, since it will produce millions of top quarks already in the first low- 
luminosity phase [1]. The great wealth of data that will be obtained at the LHC will 
allow precise measurements of the top-quark related observables. In turn, the latter must 
be matched by equally precise calculations of the relevant cross sections and differential 
distributions in perturbative QCD. 

At Tevatron, top quarks are primarily produced in pairs with their antiparticles. The 
same situation will be encountered at the LHC, where experimental collaborations antic- 
ipate measurements of the total top-quark pair-production cross section with a relative 
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error between 5% and 10%. On the theory side, the top-quark pair-production cross sec- 
tion was calculated at next-to-leading-order (NLO) in perturbative QCD [2-9]. The NLO 
electroweak corrections were obtained in [10-12]. The resummation of logarithmic terms, 
which become large near the production threshold, was extensively studied in [13-19]. 

The current resummation-improved NLO predictions for the top-quark pair-production 
cross section at the LHC shows an uncertainty of about 15% [1]; the latter is dominated 
by the scale uncertainty. Therefore, in order to reduce the theoretical uncertainty at the 
same level of the expected experimental error, the calculation of the next-to-next-to-leading 
order (NNLO) perturbative QCD corrections is required. 

The Feynman diagrams needed for the evaluation of the NNLO QCD corrections to the 
top-quark pair production can be grouped in three categories: i) two-loop corrections to the 
three-level production channels qq — > ti and gg tt, ii) one-loop matrix elements with an 
additional parton in the final state, and iii) tree-level matrix element with two additional 
partons in the final state. The last two sets of diagrams were already evaluated in the 
context of the calculation of the NLO corrections to the production of tt + Ij [20]. Contri- 
butions arising from the interference of one-loop diagrams in both the quark-antiquark and 
the gluon- fusion channels were studied in [21-24]. The first steps toward the calculation 
of the two-loop corrections were taken in [25,26]; in these papers all the relevant two-loop 
diagrams were evaluated in the limit s, |u| S> m?, where s is the squared center of mass 
energy, t is the squared momentum transfer, u = 2m?' — s — t and m is the heavy-quark 
mass. A full numerical calculation of the two-loop virtual corrections in the qq — > ti chan- 
nel was carried out in [27]. Finally, the analytic calculation of the diagrams involving a 
closed quark loop in the quark-antiquark channel was presented in [28]. 

In this paper we describe the calculation of a conspicuous subset of the two-loop planar 
diagrams in the qq — > ti production channel. The calculation of these corrections allows us 
to obtain an analytic expression for the leading color coefficient in the interference of the 
two-loop matrix element with the tree-level amplitude. Our results are valid for generic 
values of the Mandelstam invariants s, t and of the heavy-quark mass m. The calculation 
was carried out by means of a technique based on the identification of a set of Master 
Integrals (Mis) through the Laporta algorithm [29], and on their subsequent evaluation 
by means of the differential equation method [30]. The results are written in terms of a 
suitable base of one- and two-dimensional Harmonic Poly logarithms (HPLs) [31-33]. 

The paper is organized as follows. In Section || we introduce our notation and con- 
ventions. In Section ^ the method employed in the calculation is briefly described and 
the technical difficulties we met are discussed. In Section ^ we collect the explicit ex- 
pression of the counterterms needed for the ultraviolet renormalization of the Feynman 
diagrams which we evaluated. The results which we obtained and their expansion near the 
production threshold are discussed in Section Section ^ contains our conclusions. In 
the appendices we describe in some detail the HPLs we employed in the calculation and 
their expansion at production threshold. Furthermore, we collect the explicit expression of 
several previously unknown Mis belonging to box topologies. 
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Figure 1: Tree-level amplitude. Massive quarks are indicated by a thick line. 



2. Notation and Conventions 

We consider the scattering process 

q{pi) + q{P2) 



(2.1) 



in Euclidean kinematics, where pf = for i = 1,2 and pj = —m? for j = 3,4. The 
Mandelstam variables are defined as follows 



ipi + P2f , t = - {pi- psf , u = - (pi - PiY 



(2.2) 



Conservation of momentum implies that s + t + u = 2m^. 

The squared matrix element (averaged over the spin and color of the incoming quarks 
and summed over the spin of the outgoing ones), calculated in d = 4 — 2s dimensions, can 
be expanded in powers of the strong coupling constant as as follows: 



\M\'^{s,t,m,e) 



47r^a| 



\ IT 



(2.3) 



The tree-level amplitude involves a single diagram (Fig. |T]) and its contribution to Eq. 
is given by 

"(t — m^)^ + (n — m^)^ . 2m^ 



^0 = 47Vc Cf 



+ 



(2.4) 



where A^^^ is the number of colors and Cp = {N^ — l)/2Nc. 

The 0{as) term in Eq. ( |2.3D arises from the interference of one-loop diagrams with 
the tree- level amplitude [2-8]. The 0{a'g) term A2 consists of two parts, the interference 
of two-loop diagrams with the Born amplitude and the interference of one-loop diagrams 
among themselves: 



A2 = A2 



(2x0) 



■A. 



(ixl) 



The latter term ^2^^"^^ was studied extensively in [21,22]. ^2^'^^'', originating from the 
two-loop diagrams, can be decomposed according to color and flavor structures as follows: 



(2x0) 



-(2x0) 
•^2 



NcCf 



N^A + B + -^ + Ni (n,Di + ^]+Nh (n^Dh + 



' c 

+NfFi + NiNhFih + NlFh 



Nr 



Nr 



(2.5) 



-3- 




mumm 



o 
o 
o 

o. 



Figure 2: Some of the two-loop planar box diagrams involved in the calculation. 



where A'^; and Nh are the number of Hght- and heavy-quark flavors, respectively. The coef- 
ficients A, B , . . . , Fh in Eq. ( |2.5| ) are functions of s, t, and m, as well as of the dimensional 
regulator e. These quantities were calculated in [25] in the approximation s, |u| m^. 
For a fully differential description of top quark pair production at NNLO, the complete 
mass dependence of ^2^^^^^ is required. An exact numerical expression for it has been 
obtained in [27]. Exact analytic formulae for all the coefficients arising from Feynman 
diagrams involving closed quark loops {Di, Ei, D^, Eh,Fi, F^, Fi^) were obtained in [28]. In 
this work, we provide an exact analytic expression for the coefficient A in Eq. (^]^), which 
arises from planar Feynman diagrams only. 



3. Calculation 

The package QGRAF [34] generates 218 two-loop Feynman diagrams contributing to the 
process qq tt. After treating the QGRAF output with FORM [35] in order to carry 
out the color (and Dirac) algebra, one finds that there are 44 non-vanishing diagrams con- 
tributing to the leading color coefficient in Eq. (|2.5| ). Some of the box diagrams involved 
in the calculation are shown in Fig. ^. All the two-loop graphs encountered in our calcula- 
tion can be treated by employing the technique based upon the Laporta algorithm for the 
reduction to a set of Master Integrals (Mis) . The Mis are then evaluated by means of the 
differential equation method. 

A subset of the Mis needed in the calculation were available in the literature [36-41]. 
They were employed in previous two-loop calculations of the heavy-quark form factors [42] , 
amplitudes for Bhabha scattering [43], heavy-to-light quark transitions [44], and in the 
calculation of the fermionic corrections to qq tt [28]. The Mis not included in this 
subset are a part of the original findings of this paper. 

The reduction to Mis was carried out with an implementation of the Laporta algorithm 
written in C++ [45], and large parts of it were cross checked with the Maple package A. I .R. 
[46]. The six so far unknown irreducible topologies encountered in the calculation of the 
planar diagrams are shown in Fig. ^. In the figure, thick internal lines indicate massive 
propagators, while thin lines indicate massless ones. An external dashed leg carries a 
squared momentum (pi + P2)'^ = —s; other external lines indicate particles on their mass- 
shell, where pf = for thin lines and pf = —m? for thick lines. The Mis in Fig. |^ are 
collected in Appendix ^ and in a file included with the arXiv submission of the present 
paper. 
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(a) - 2 Mis 



(b) - 2 Mis 



(c) - 2 Mis 




(d) - 2 Mis 




(e) - 2 Mis 




(f) - 3 Mis 



Figure 3: New non-reducible topologies encountered in the calculation of the planar diagrams. The 
number of Master Integrals related to each topology is indicated in the figure. 



In calculating the Mis by means of the differential equation method, it is crucial to 
fix the undetermined integration constant(s) through a boundary condition. While there 
is no general method available to fix this boundary condition, it is usually sufficient to 
know the behavior of the MI at some particular kinematic point. For example, knowing 
that the integral is regular for a certain value of s, one can impose the regularity of the 
solution of the differential equation at that point. In cases in which considerations related 
to the regularity of the solutions of the differential equations are not sufficient to fix all of 
the integration constants, it is necessary to directly evaluate the Mis at a specific phase- 
space point. This can be done using techniques based on Mellin-Barnes representations 
of the integrals [47], with the help of the Mathematica packages Ambre [48] and MB [49]. 
In order to numerically check the analytic calculation of the Mis, we employed the sector 
decomposition technique [50], implemented in the Mathematica package FIESTA [51]. 

All the Mis were calculated in the non-physical region s < 0, where they are real and 
can be conveniently written as functions of the dimensionless variables 

+ Am? - -t ^^^^ 



V-s + 4m2 + ' m? ' 

The transcendental functions appearing in the results are one- and two-dimensional har- 
monic polylogarithms (HPLs) [31-33, 40] of maximum weight four. In Appendix ^ we 
briefly review the definition of the HPLs appearing in the calculation. All the HPLs ap- 
pearing in the analytic expression of the coefficient A can be evaluated numerically with 
arbitrary precision by employing the methods and codes described in [32]. 

Following the procedure outlined in the present section, it is possible to obtain the 
expression of the bare squared matrix elements involving planar two-loop diagrams. The 
renormalization of the ultraviolet divergencies is discussed in the next section. 



4. Renormalization 



The renormalized QCD matrix element is obtained from the bare one by expanding the 
fohowing expression : 

•^ren = -^wf,?! ^^bare ("s.bare ^as^^S , ""T-bare ~^ ^mm^ , (4.1) 

n 

where ZwY,n is the external leg wave function renormalization factor, as is the renormalized 
coupling constant and m is the renormalized heavy-quark mass. (In the rest of the section 
we suppress the subscript "5" in as)- 

We postpone the discussion of mass renormalization to the end of the section and we 
start by considering the coupling constant and wave function renormalization. 

We introduce the following quantities: 

ao = ^^, and a = -. (4.2) 
vr TT 

By expanding the amplitude and the wave function renormalization factor in oq we find: 

Aen(abare) = "oA + a^Ai + + ©(oq) , 

ZwF,„ = 1 + aoSZ^l^ + aldZ^l^ + Oia^) . (4.3) 

The relation between oq and a is given by: 

ao = a + a^6z2^ + a^6Zj^^ + O(a^) . (4.4) 

By employing Eqs. (4.3,4.4) in Eq. (4.1) we find 

Aen = aAo + a^^^-en + oMj-en + C(a^) , 

4eli = -^i+{Y. l^Z^ln + SZ2^^ Ao , 

+ \sz!^l,n + SZi'^ Yl '^^wL + <54'^ ) ^0 • (4.5) 

n n / 

In the equations above, Ai represents the bare amplitude at i loops stripped of the fac- 
tor a. In the case of the process qq — > ti, the wave function renormalization factors of 
massless quarks vanish at one loop, while the ones of the massive quarks in the on-shell 
renormalization scheme are given by 



my ^\ 4.e l-2e 



5Z'^' = C{e) {^] Cp{--- ) , (4.6) 



where the subscript M indicates massive quarks and where C{e) = (47r)'^r(l + e). The 
one-loop renormalization constant for a in the MS scheme is given by 
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Figure 4: One-loop diagrams (excluding the diagrams with closed quark loops). Thin arrow lines 
represent massless quarks, thick arrow line massive quarks, dashed arrow lines are Faddeev-Popov 
ghosts, and coiled lines are gluons. 



where Pq = 11/6Ca — 1/3(A^/ + Nh) and 7 is the Euler-Mascheroni constant 7 
Therefore, the overall one-loop counter term is given by 



C{s) 
4e 



2po + 3 + 4e + 3eln 



Cf + O (e^) 



0.577216. 



(4.8) 



Here we are interested in finding the ultraviolet counterterm needed to renormalize the 
two-loop diagrams not involving closed quark loops. The part of this counterterm needed 



to renormalize the leading color structure in Eq. (2^) can then be trivially extracted. 
Taking into account the fact that the wave function renormalization factors are zero for 
the incoming particles and identical for the massive ones, one finds that the renormalized 
amplitude (excluding quark-loop diagrams) is 



CT ) 



1=5 



WF,M 



(4.9) 



In Eq. ( |4.9| ), the quantity A^i^^ is the amplitude of the j-th diagram in Fig. ^ (stripped 
of the factor a). The quantity " ) indicates the ^-th diagram in Fig. ^ with a 

mass counter term insertion in one of the internal heavy quark lines. The renormalization 
coefficients not previously defined are: 



-Cie) 
C{ef 



r(l + e) ^12e 



-7e 



r(l + e) 



1 

4e 



r2 



17 

— ( 

12 
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9 51 433 3,, , n, 13 o 

+ + C 3) + vr^ In 2 vr^ 

32^2 64e 128 2^^ ^ 16 . 



^ , 11 101 803 3^, , vr^ 5 o > . x 

+Ca ^ + -C 3 n 2 + — TT^ + O e 

* 32e2 64e 128 4^^ ^ 2 16 ' 



(4.10) 



and they can be found in [25,52]. 
5. Results 

The main result of the present paper is an analytic, non-approximated expression for the 
coefficient A in Eq. ( ^.5] ). Since such a result is too long to be explicitly printed here, we 
included in the arXiv submission of this work a text file with the complete result, which 
is written in terms of one- and two-dimensional HPLs of maximum weight four. Since the 
coefficients in Eq. ( |2.5D still contain infrared poles, the result is dependent on the choice 
of a global, e-dependent normalization factor. With our choice, we factor out an overall 
coefficient 

C\e) = [{4nrril + e)]\ (5.1) 

We also provide a code that numerically evaluates the analytic expression of the quantities 
listed above for arbitrary values of the mass scales involved in the calculation. The code is 
written in C++ and uses the package for the evaluation of multiple polylogarithms within 
GiNaC [32]. 

In order to cross check our results, we expanded them in the s, \t\, \u\ ^ m? limit. The 
first term in the expansion agrees with the results published in [25] ; the second order term 
agrees with the results found in the Mathematica files included in the arXiv version of [27]. 
We also find complete agreement with the numerical result of Table 3 in [27], corresponding 
to a phase-space point in which the s, \t\, \u\ ^ m? approximation cannot be applied. With 
our code it is also possible to reproduce the first plot in Figure 4 of [27], where the finite 
part of A is shown as a surface depending on the variables rj and defined as 

* = -^. Ul-.^<^<Ul^.^. (5.2) 



4m2 ' ^ s ' 2V y l + r]J - ^ - 2 \ ]/ I + r] 

This surface is shown on the left side of Figure ^. 

Furthermore, it is possible to expand our result for values of the center of mass energy 
close to the production threshold. We define 



.„.,.(^). ,.3) 

where 9 is the scattering angle in the partonic center of mass frame, and we expand our 
results in powers of the heavy quark velocity (3, up to terms of order f3'^. The coeffi- 
cients of this expansion contain transcendental constants which originate from one- and 
two-dimensional HPLs evaluated at x = 1. Since we did not find a satisfactory analyti- 
cal representation for all of these constants, in the formulae below we present them in a 
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Figure 5: Left: finite part of the coefficient A as a function of the variables rj and (j). Right: the 
exact value of A^^^ as a function of j3 (red solid curve) versus its expansion close to threshold up to 
terms of order (blue dashed curve). Both curves are plotted for ^ = 1/2. In both cases we used 
the normalization adopted in [27] to facilitate comparisons. 



numerical form. We find: 

A(-^\^,o A(-^H^,0 a^'^Hp,^ 



^(-4) 
^(-2) 



+ 



+ 



+ 



0.25 - 0.5/3^ (1 - e + O , 

1.68185 + 0.5L^ + /3(1 - 2^ - [o.5 + ^(1 - 0(5.86371 + L^)] + O {(3^) , 



-2.67119 - 0.302961L„ + O.bLl + (3 



1.84475 + 2L^ - ^3.68951 + 4L^ 



0.777936 -Lf,- ^(1 - 0(7.03784 + 4.39408L^ + L^) +0 {(3^) 



yl(-i) = -8.15701 - 5.7593L^ - 2.13629L2 + 0.333333L^^ + /3 



-9.83935 - 3.64382L, 



+2Ll + ^(19.6787 + 7.28765L^ - 4L^ 



2.78693 + 5.22254L, 



-Ll + ^(1 - 0(46.7006 + 8.75816L^ - 0.727411L2 - 0.666667L^)j + O (/3^) , 



23.5701 + 7.82592L^ + 0.754463L^ '^'"^''^ 



2.03531Lf, + 0.166667L2 + 



0.505501 - 11.5953L^ - 7.31049L2 + 1.33333L|^ ^-l-OH + 23.1906L^ 



+14.621L2 - 2.66667L|^) + 



-4.351 - 3.60348L„ + 7.05587L2 



-0.666667L)^ + ^{1 - 0(-5.36823 + 43.1864L^ + 6.73063^^ + 0.737281L)^ 



-0.333333L^] 



(5.4) 



Note that the dependence on (3 and on ^ in the formulae above is only polynomial. All 
the logarithmic terms ln/3, In^, ln(l — 0) lii(l ~ • • •) which are indeed present in the 
expansion of individual HPLs, cancel out in the final expressions. The coefficient A is finite 
at threshold. The expansion presented here could be used in the future for the calculation 
of logarithmically enhanced terms near the ti production threshold. On the right hand 
side of Fig. ^, we compare the exact expression of the coefficient A^'^'f with the expansion 
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in powers of /? including up to terms of order 0^ (in the plot we set ^ = 1/2). 
6. Conclusions and Outlook 

In this paper, we presented the analytic calculation of the two- loop planar corrections to 
the heavy-quark production amplitude for qq tt, retaining the exact heavy-quark mass 
dependence. The calculation of this subset of diagrams allowed us to obtain an analytic 
expression for the leading color coefficient in Eq. ( |2.5D . We checked our formula against 
recent results obtained in analytic form in the small-mass region [25]. Moreover, we found 
numerical agreement with the value of A at the phase-space point presented in [27] . 

Our result represents a gauge invariant sub-set of the full two-loop corrections to 
the partonic process qq — > ti. In order to complete the analytic calculation of the two- 
loop corrections, it is necessary to calculate the non-planar diagrams. Likewise, analytic 
results for the two-loop amplitude for gg — > tt could be obtained in the same calculational 
framework [53]. However, some of the two- loop diagrams appearing in the gluon fusion 
channel cannot be expressed in terms of two-dimensional HPLs. In fact, their reduction 
to Mis involves a "sunrise" -type subtopology with three equal massive propagators and an 
external momentum which is not on the mass shell of the internal propagators. It is known 
that already such a three-propagator graph involves elliptic integrals [54]. 

In order to obtain NNLO predictions for the total ti production cross section and 
for differential distributions, it is necessary to combine the two-loop virtual corrections 
with the already available [20] one- loop corrections to the tt+{l parton) process and with 
the tree-level tt+{2 partons) process. These diagrams with additional partons in the final 
state contribute to infrared-divergent configurations where up to two partons can become 
unresolved. Their implementation requires the application of a NNLO subtraction method. 
The methods presently available [50, 55, 56] have been applied up to now [57-60] to at 
most 1 — > 3 processes in e^e~ annihilation and 2 — > 1 processes at hadron colliders. A 
calculation of a hadronic 2 — > 2 process, involving massive partons, will represent a new 
step in complexity, potentially requiring further refinements of the methods available to 
date. 
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A. Harmonic Polylogarithms 

The results presented here are conveniently expressed in terms of one- and two-dimensional 
HPLs. Nowadays, harmonic polylogarithms are extensively employed in multiloop com- 
putations. Therefore, in this appendix we only briefly review their definition. The reader 
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interested in the algebraic properties of these functions can find detailed discussions about 
this topic in the available literature [31-33,40]. 

In the non-physical region s < 0, seven weight functions are needed for the HPLs with 
argument x. They are^ 

fu,{x) = , with u;e <^ 0,1,-1, -y,—,-±^ ^ . A.l 

X — w I y 1 1 \ 

For HPLs with argument y, we need six weight functions 

fw{y) = —^, with |o,l,-l,-x,--,l---xl . (A.2) 
y — w I XX) 

The weight-one HPLs are defined as 

G(0;x)=lnx, G{w;x) = [ dtf^{t) . (A.3) 
HPLs of higher weight are defined by iterated integrations 

G{w,----x)= r dtU{t)G{---;t), (A.4) 



Jo 

with the only exception being the HPLs in which all the weights are zero which are defined 
as follows 

(7( 0, 0, y ,0 ; x) = — In" a; . (A.5) 

n 

Analogous definitions hold for HPLs with argument y. The reader should be aware of 
the fact that in the original definition of Remiddi and Vermaseren, the weight function 
corresponding to the weight 1 was /i = 1/(1 — x). In order to translate the HPLs defined 
with the Remiddi- Vermaseren convention to the ones employed in this work (and vice 
versa) it is sufficient to multiply each HPL by a factor (—1)"", where n is the number of 
weights equal to 1. 

The weights —y and —1/y were already introduced in [28,33,43]. In our results, the 
two-dimensional harmonic polylogarithms have maximum weight four. Therefore, it was 
not possible to rewrite all the two-dimensional HPLs in terms of Nielsen polylogarithms, 
as it was done for the results of [28] , where the two-dimensional HPLs had maximal weight 
three. However, it is possible to evaluate all the HPLs appearing in the analytic expres- 
sion of the coefficient A in Eq. (|2.5| ) by employing the GiNaC implementation of multiple 
polylogarithms by Vollinga and Weinzierl [32]. 

We first obtained the squared matrix elements in the non-physical region s < 0. The 
corresponding quantities in the physical region s > 4m^ could be obtained by analytic 



The last two weie hts in Eq. ( |a1| ) introduce explicit imaginary parts in the formulae. However, these 
HPLs appear in such a way that these imaginary parts cancel in the non-physical region, where the result 
has to be real. Alternatively, one could choose a basis of weight functions in which the pair i ± is 
replaced by the original quadratic expressions in the integrals: l/{x^ — a; + 1) and x/{x^ — x + 1) [40]. In 
this case these HPLs are all manifestly real. 
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continuation to the complex value s — > s + i(5, where 6^0'^. For s > Am? the variable x 
becomes 

x = -x' + i5, (A.6) 

where 

x' = , (A.7) 

y/s + V s — Am? 

So that < x' < 1 for 4m^ < s < oo. The HPLs of argument x can develop an imaginary 
part because of analytic continuation^. In particular, the imaginary part of the HPLs of 
argument x for s > Am^ is defined when the analytic continuation of the logarithm is 
specified: 

G(0; x) = G(0; -x' + i5) = G(0, x') + m . (A.8) 
For notation convenience, after the analytic continuation we rename x' as x. 

B. Expansion of the HPLs near the Threshold 

We devote this section to a brief discussion of the technique employed to expand the 
coefficient A near the production threshold (3 = ^Jl — Am^ / s = 0. 

The first step consists in carrying out the analytic continuation from the non-physical 
region, s < 0, to the physical one, s > 4m^, according to the method outlined above. The 
one- and two-dimensional HPLs appearing in the analytically continued expression of the 
coefficient A must then be expanded in the /? — > limit. While the threshold expansion 
of the ordinary HPLs does not lead to any particular difficulty, the expansion of the two- 
dimensional HPLs is indeed more delicate. The reason is that, in the latter case, the 
expansion parameter f3 appears in both the argument and the weights of the HPLs. In 
fact, in the physical region one finds: 

Moreover, the coefficient A depends on two-dimensional HPLs of maximal weight four; 
therefore it is very challenging to obtain explicit analytic expressions of these functions in 
terms of logarithms and poly logarithms of complicated arguments, which can be subse- 
quently expanded in /? — > 0. Such an approach should be replaced by a more direct and 
algorithmic method. 

In the following we describe the method which allows to extract directly the expansion 
of a given two-dimensional HPL of weight n, assuming that the expansion of the HPLs 
of weight n — 1 is known. Let us consider, for simplicity, the case in which the HPL has 
argument x and it has y (or 1 /y) in the weights. Since the dependence of x on /3 is the one 
shown in Eq. ( |B.1| ), we first use the transformation that relates G{w, (1 — /3)/(l -|- to 
G{w' , (3). This transformation allows to rewrite the HPL function of x as a combination 



^The coefficient A is real for s < 0. However, because of the weight functions we use, and because of the 
fact that < J/ < oo, individual HPLs can develop imaginary parts also in the non- physical region. The 
latters cancel out among each other. 
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of HPL functions of /3, HPLs with y in the weights but evaluated in x = 1, and HPLs 
(either one- or two-dimensional) of a smaller weight. The series expansion of the HPLs 
with argument P is found recursively. We write them as the integral between and P of 
the total derivative with respect to /3 of the HPLs themselves. The total derivative gives 
rise to HPLs of lower weight. We insert the expansions of the lower- weight HPLs and then 
we integrate again. For the HPLs evaluated in j; = 1 the procedure is analogous. It can 
happen that in the intermediate steps (expansion followed by an integration) logarithmic 
divergences occur. These divergences must be regularized and they cancel in the final 
expressions. 

Let us illustrate the algorithm more in detail. To this purpose, we consider the example 
of a simple two-dimensional HPL of weight two, which appears in the expression of the 
coefficient A. For < x < 1 and x < y < 1/x we define 

G(y,-l;x)= dt G(-l,t) = dt ln(l + 1) . (B.2) 

Jo t-y Jo t-y 

This HPL is real in the physical region, however the intermediate stages of the procedure 
described below require the sum of complex terms. We assume that any ambiguity of this 
sort is dealt with by assigning an infinitesimal imaginary part to y. We start by rewriting 
the HPL of argument x in terms of HPLs of argument (3; the integral representation in 
Eq. ( |B.2| ) can be rewritten as 

G{y, -l;x) = G{y, -1; 1) + T dt-^G{-l, t) . (B.3) 

J I t-y 

In the integral in the second term of Eq. ( |B.3| ) we carry out a change of integration variable 
by setting q = {I - t)/{l + t) {dt = -2/(1 + q)'^dq) 

'dt^Gi-l,t) ^ -2 rdq^J^±4^'G(-l,' 



1 t-y ' Jo ^(l-g)-y(l + g)(l+g)2 ^ 'l + 



/3 







1 1 \ ^ / . 1 



G -1;--^ , (B.4) 



-7 g + iy V 1 + 

where we employed the relation f3 = {1 — x)/{l + x) and we introduced 7 = (1 — y)/(l + y). 
We emphasize that, given the definitions of x and y in the physical region, 7 is a function 
of P and of the variable ^ defined in Eq. ( ^.Sj) . The weight one HPLs appearing in the 
integrand of Eq. (B.4) can also be easily rewritten in terms of HPLs of argument q: 



'^(-I'lT^) =-G(-l;g)+ln2. (B.5) 
At this stage it is straightforward to integrate in q to finally obtain 

Giy, -1; x) = G{y, -1; 1) - (^(7, -1; /?) - G(-l, -1; /?)) 

+ ln2fG(7;/3)-G(-l;/3)V (B.6) 



In the equation above, there are one- and two-dimensional HPLs of weight one, that we 
consider to be known, and three HPLs of weight two. Among them, there is a single two- 
dimensional HPL of weight two, G —1; P) , and two one-dimensional HPLs. One of 



- 13 - 



them, G(— 1, — 1; /3), has a trivial expansion in /3 — > 0. Let us discuss the method employed 
to obtain the threshold expansion of the other two: G (7(/3), — 1; /?), and G{y, —1; 1). 
We can rewrite G (7(/3), — 1; /3) as follows: 



G (7(/3), -1; /?) = dp'^G (7(/3'), -1; P') + G (7(0), -1; 0) . 



(B.7) 



In this simple example the second term in Eq. ( |B.7| ) is well defined (and it is actually equal 
to zero) . The derivative in the first term of Eq. ( |B.7| ) can be rewritten as 



d ^ 
dp* 



13 


-7(/3) 




1 




-7(/3) 




1 


P 


-liP) 




1 


P 


-liP) 




1 



l + 7(/3) 



G{-l;P)- 
G{-l;P)- 
dt 



dt 



dp h (t-iiPW^ ' ' 

dp r^4(^)^(-^'^) 



d-f{P) ff^ , d 



drm 
dp 

1 



1 



P-l{P) 
1 



G{-l;P) 



t - 7(/3) t + 1 



-G(-l;/3) 



/3-7 



G(-l;/3) 



1 + 7 



G(7;/3) 



(B. 



where in the last line we dropped the dependence of 7 on p. Since the expansion of the 
HPLs of weight one is assumed to be known, it is straightforward to expand the equation 
above in the limit /? — > and to insert it in Eq. ( |B.7] ) to obtain 



G(7,-l;/?) =/? l + (l-20 ln2 + ln(0-ln(2e-l) +0{p' 



(B. 



The formula above is real for ^ > 1/2. However, the imaginary parts which arise for ^ < 1/2 
cancel against the imaginary parts coming from the expansion of G{y, —1; 1). 

The expansion of G{y, —1; 1) can be done with the same algorithm. We write: 

G{yiP), -1; 1) = r 4^,G{y{P'), -1; 1) + G(y(0+), -1; 1) , (B.IO) 
7o+ dp' 

where 0^ indicates the fact that we must take the limit /? ^ 0"^ both in the integration 
constant and in the lower boundary of the integration. Both limits are logarithmically 
divergent, but the divergencies cancel between the two terms. We have: 



dt 



dyjp) 
dp Jo {t-y) 



rG(-l,t) 
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dp 
dp 



1 



t-y 

In 2 



G{-l,t) 
1 

+ 



dt 



t-y 1 + 



1-2/ 1 + y 



G(y;l)-ln2 



(B.ll) 



Expanding Eq. ( [B.ll| ) in the limit /3 — > 0, we find the following Laurent series: 
^G{y{P'), -1; 1) = ^ + ^ + (2e - 1) [In + In (2^ - 1)] + 0{P) . (B.12) 
Integrating Eq. ( B.12 ), as in Eq. ( |B.10| ), we find: 



^ TRl^^y^^'^^ -1; 1) = In 2 In - In 2 In (0+) + pl,^ + (2^ - 1) In 



+ ln(2e-l)-l \ + 0{p^). 



(B.13) 



Moreover, one finds: 



G(y(0+), -1; 1) = G(y(0+); 1)G(-1; 1) - G(-l, y(0+); 1) , 

= In 2 [In (0+) + In 2 + In (2^ - 1)] + \c,2-\ In^ 2 . (B.14) 



Finally, we have: 

G{y{P), -1; 1) = In 2 In (/?) + In 2 In (2^ - 1) + 1 In^ 2 + 



-P 



In 2 
2^ - 1 



+ (2e - 1) In (/?) + In (2e - 1) - 1 + ©(Z^^) , (B.15) 



where the divergencies disappeared. Considering also the expansion of the one- and two- 
dimensional HPLs of weight 1 and the expansion of G(— 1, —1; P) in Eq. ( |B.6| ), we find the 
final formula: 

= ln21n(/3)+ln21n(e) + ^ln2 2 + ^C2 

+p[-2i + (2e - l)[ln (/?) + In + ln2]} + 0{P^) . (B.16) 

In the case in which the argument of the HPLs is y and x is present in the weights, 
G{w, the procedure is analogous to the one explained above. Since the dependence 



of y on /? involves also the parameter ^ (see Eq. (B.l)), the first step consists in using the 
scale properties of the HPLs: 



G{w,...;y) = G{Xw,...;Xy) , 



(B.17) 



(valid in the case in which all the trailing zeroes have already been extracted) to get rid 
of the ^ dependence in the argument. To achieve this goal, one multiplies weights and 
argument by 



A 



1 + 



4/3 
(1-/3) 



re 



(B.18) 
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In so doing, we fit again in the case illustrated in the example concerning G{w, x), since 



and the expansion proceeds along the same steps as outlined above. 
C. Master Integrals 

In this Appendix we collect the Mis for the topologies in Fig. ^. 

The explicit expression of the Mis depends on the chosen normalization of the inte- 
gration measure. The integration on the loop momenta is normalized as follows: 



where C(e) was defined in Eq. In Eq. ( |C.1| ) fi stands for the 't Hooft mass of 

dimensional regularization. The integration measure in Eq. ( p.l| ) is chosen in such a way 
that the one-loop massive tadpole becomes 

2^ ^ 72 2=-7n — r- C.2 

^2 _^ ^2 _ 

In calculating the squared matrix element, we multiply our bare results by (/i^/m^)^, in 
order to make explicit the dependence on the 't Hooft scale. We also point out that, since 
the squared matrix element still contains soft and collinear divergencies regulated by e, it 
depends on the normalization of the integration measure. In particular, in order to match 
our results with the ones of [25,27], it is necessary to multiply the latter by the factor 

= 1 _ ^(2)^2 + 2^(3),3 + ^Ci2?e' + O {e') . (C.3) 

The Mis are expanded in powers of the dimensional regulator e; below we collect the 
analytic expression of the coefficients in the e expansion up to terms involving HPLs and 
2dIIPLs of weight three. The coefficients involving HPLs and 2dHPLs of weight four 
are also needed in order to obtain the finite part of the leading color coefficient, and 
we calculated them. However, their analytic expressions are too long to be written in 
this appendix; the interested reader can find them in the text file included in the arXiv 
submission of this paper. 

There are two Mis belonging to the topology Fig. ^(a). The first MI is 




^"fel^"fc2 

PO {k2) Po ih - k2) Po {k2 - Pi) Po {k2 -Pi- P2) Pm {kl " Ps) ' 

where we define 

Po (k) = k\ Pm (k) = k^ + m'^. (C.5) 
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We find 




(C.6) 



^-3 
A-2 



X 



16(1 -x)^ 

X 



A. 



16(1 — x)'^y 

X 



1 



-2y - yG{0; x) + 2yG(l; x) + {y + 1)G(-1; y) 
y(3C(2) - 4) + 2{y + 1)G(-1; y) - 2yG(0; x) 



16(1 — x)^y 

+{y + 1)G(-1; y)G(0; x) + 4yG(l; x) - 2(y + 1)G(-1; y)G{l; x) 
-2{y + 1)G(-1, -l-y) + {y- 1)G(0, -1; y) - yG(0, 0; x) + 2yG(0, 1; x) 
+22/G(l,0;x) -4yG(l,l;x) 



A. 



X 



2y(3C(2) + 7C(3) - 4) - (7C(2)y - 4y + 11C(2) - 4)G(-1; y) 



16(1 -x)2y 

-y(C(2) + 4)G(0; x) + 2{y + 1)G(-1; y)G(0; x) - 2y(3C(2) - 4)G(1; x) 

-4(y + 1)G(-1; y)G(l; x) + 4(y + 1)C(2)G (-1/y; x) - 4(y + 1)G(-1, -1; y) 

-4(y + 1)G(0; x)G(-l, -1; y) + 4(y + 1)G(1; x)G(-l, -1; y) 

+2(y + 1)G (-1/y; x) G(-l, -1; y) + 2(y + l)G(-y; x)G(-l, -1; y) 

+2(y - 1)G(0, -1; y) + 2yG(0; x)G(0, -1; y) - 2(y - 1)G(1; x)G(0, -1; y) 

-{y + 1)G (-1/y; x) G(0, -1; y) - (y + l)G(-y; x)G(0, -1; y) - 2yG(0, 0; x) 

+4yG(0, 1; x) + 4yG(l, 0; x) - 2(y + 1)G(-1; y)G(l, 0; x) - 8yG(l, 1; x) 

+4(y + l)G(-l;y)G(l,l;x) + (y + l)G(-l;y)G(-l/y,0;x) 

-2(y + l)G(-l;y)G(-l/y,l;x) + (y + l)G(-l;y)G(-y,0;x) 

-2(y + l)G(-l;y)G(-y,l;x) + 2(y-l)G(-l,0,-l;y)+4G(0,-l,-l;y) 

-2yG(0, 0, -1; y) - 2yG(0, 0, 0; x) + 4yG(0, 0, 1; x) + (3y + 1)G(0, 1,0; x) 

-2(3y + 1)G(0, 1, 1; x) + 2yG(l, 0, 0; x) - 4yG(l, 0, 1; x) - 4yG(l, 1,0; x) 

+8yG(l, 1, 1; x) + (y + 1)G (-1/y, 0, 0, x) - 2(y + 1)G (-1/y, 0, 1; x) 

-(y + 1)G (-1/y, 1, 0; x) + 2(y + 1)G (-1/y, 1, 1, x) - (y + l)G(-y, 1, 0; x) 

+2(y + l)G(-y,l,l;x)] . (C.7) 

Tlie second MI for the topology in Fig. |3|-(a) was chosen as 




Po {k2) Po ih - k2) Po {k2 - Pi) Po {k2 -Pi- P2) {h " Ps) 



(C.8) 



Its analytic expression is 
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1 ° 



(C.9) 



A. 



A. 



G'(-l;y), 



An 



y)G(0; x) + 2G(-1; y)G(l; x) + 2G(-1, -1; y) 



11C(2)G(-I;y) + 2G(1, 0; y) - 4G(1, l;x)G(-l;y) 



16(1 — 

X 

16(1 — x)^y 
+G(0,-l;y)" 

X 

16(1 — x)^y 

-G (-l/y, 0; x) G(-l; y) + 2G (-l/y, 1; x) G(-l; y) - G{-y, 0; x)G(-l; y) 
+2G(-y, 1; x)G(-l; y) - 4C(2) G (-l/y; x) + 4G(0; x)G(-l, -1; y) 
-4G(1; x)G(-l, -1; y) - 2G (-l/y; x) G(-l, -1; y) - 2G(-y; x)G(-l, -1; y) 
-2G(1; x)G(0, -1; y) + G (-l/y; x) G(0, -1; y) + G(-y; x)G(0, -1; y) 
+2G(-1, 0, -1; y) - 4G(0, -1, -1; y) - G(0, 1, 0; x) + 2G(0, 1, 1; x) 
-G(-l/y,0,0;x)+2G(-l/y,0,l;x) + G(-l/y,l,0;x) 



-2G (-l/y, 1, 1; x) + G(-y, 1, 0; x) - 2G(-y, 1, 1; x) 
The topology shown in Fig. ^(b) has two Mis. The first one, 



(C.IO) 




PO ih) Po {h - k2) Po {k2 - Pi) Po (kl -PI- P2) Pm {kl - P3) 



, (C.ll) 



has the following analytic expression: 




(C.12) 



i=-l 



X 



[4C(2) + G(0,0;x)-2G(0,l;x) 



16(1 - x2 

^^^^"1^2) [C(2) (l6G(-l; x) - G(0; x) - 12G(1; x) - 4 G (-l/y; x) + 8 

+7C(3) - 2 G (-l/y; x) G(-l, -1; y) + 2G(-y; x)G(-l, -1; y) 

+G(0; x)G(0, -1; y) + G (-l/y; x) G(0, -1; y) - G(-y; x)G(0, -1; y) + 2 G(0, 0; x) 

-4G(0, 1; x) - G(-l; y)G (-l/y, 0; x) + 2G(-1; y) G (-l/y, 1; x) 

+G(-1; y)G(-y, 0; x) - 2G(-1; y) G(-y, 1; x) + 4G(-1, 0, 0; x) - 8G(-1, 0, 1; x) 

+2G(0, -1, -1; y) - G(0, 0, -1; y) + 2 G(0, 0, 0; x) - 4G(0, 0, 1; x) - 4G(0, 1, 0; x) 
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+8G(0, 1, 1; x) - 3G(1, 0, 0; x) + 6 G(l, 0, l;x) -G {-1/y, 0, 0; x) 

+2G (-1/y, 0, l;x)+G (-1/y, 1, 0; x) - 2G (-1/y, 1, 1; x) - G{-y, 1, 0; x) 

+2G(-y,l,l;x)] . (C.13) 

The second MI for the topology in Fig. |3|-(b) is: 




PO ih) Pi {h - h) Po {k2 - Pl) Po {ki -Pi- P2) Pm ih - P3) 



, (C.14) 



where 




1 ° 
— ^A,e^ + 0(e) 



(C.15) 



i=-3 



^-3 
A.2 

A-i 
A. 



7x 



192(1 -x)2(l + y) 



96(1 -x)2(l + y) L 

X 



24(1 -x)2(l +y) L 
+2G(-l,-l;y) 

X 



-8G(-1; y) + 3G(0; x) - 6G(1; x) 



-5C(2) - 3G(-1; y)G(0; x) + 6G(-1; y)G(l; x) 



-29C(3) + 46C(2)G(-1; y) + 3C(2)G(0; x) 



48(1 -x)2(l +y) 
+30C(2)G(1; x) - 36C(2)G (-1/y; x) + 24G(0; x)G(-l, -1; y) 
-12G(l;x)G(-l,-l;y) -18G(-l/y;x)G(-l,-l;y)-18G(-y;x)G(-l,-l;y) 
-9G(0; x)G(0, -1; y) + 9G (-1/y; x) G(0, -1; y) + 9G(-y; x)G(0, -1; y) 
+18G(-1; y)G(l, 0; x) - 36G(-1; y)G(l, 1; x) - 9G(-1; y)G (-1/y, 0; x) 
+18G(-1; y)G (-1/y, 1; x) - 9G(-1; y)G(-y, 0; x) + 18G(-1; y)G(-y, 1; x) 
+32G(-1, -1, -1; y) - 18G(-1, 0, -1; y) - 18G(0, -1, -1; y) + 9G(0, 0, -1; y) 
+9G(1, 0, 0; x) - 18G(1, 0, 1; x) - 18G(1, 1, 0; x) + 36G(1, 1, 1; x) 
-9G (-1/y, 0, 0; x) + 18G (-1/y, 0, 1; x) + 9G (-1/y, 1, 0; x) 

-18G (-1/y, 1, 1; x) + 9G(-y, 1, 0; x) - 18G(-y, 1, 1; x)] . (C.16) 
The are two Mis for the topology in Fig. |3|-(c). The first one is 




Po ik2) Po (kl - k2) Po {k2 - Pl) Po ih -Pl- P2) Pm ih - Ps) 



, (C.17) 



and the corresponding Laurent expansion in e is 
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A. 



C(3) + C(2)G(-l;2/)-3C(2)G(0;x) 



16(1 — X + + xy) 
-6C(2)G(1; x) + 4C(2)G (-1/y; x) - 4G(1; x)G(-l, -1; y) 

+2G {-l/y- x) G(-l, -1; y) + 2G(-y; x)G(-l, -1; y) + G(0; x)G(0, -1; y) 

-G (-l/y; x) G(0, -1; y) - G(-y; x)G(0, -1; y) - 2G(-1; y)G(l, 0; x) 

+4G(-1; y)G(l, 1; x) + G(-l; y)G (-l/y, 0; x) - 2G(-1; y)G (-l/y, 1; x) 

+G(-1; y)G(-y, 0; x) - 2G(-1; y)G(-y, 1; x) + G(-l, 0, -1; y) 

+2G(0, -1, -1; y) - G(0, 0, -1; y) - G(0, 0, 0; x) + 2G(0, 0, 1; x) 

-G(l,0,0;x) +2G(1,0, l;x) +2G(l,l,0;x) -4G(l,l,l;x) 

+G (-l/y, 0, 0; x) - 2G (-l/y, 0, 1; x) - G (-l/y, 1, 0; x) 

+2G (-l/y, 1, 1; x) - G(-y, 1, 0; x) + 2G(-y, 1, 1; x 



(C.19) 



The second MI for the topology in Fig. |3|-(c) was chosen to be 




Po ik2) Po (ki - k2) Po {k2 - Pi) Po {ki - pi - P2) Pi {ki - p^) 



, (C.20) 



where 




(C.21) 



i=-3 



^-3 
A-2 

A-i 
A_. 



X 



96(1 -x)2 

X 

96(1 - x)2 

X 



2G(-1; y) - 3G(0; x) + 6G(1; x) 

11C(2) - 6G(-1; y)G(0; x) + 12G(-1; y)G(l; x) - 4G(-1, -1; y) 



96(1 - x)2 

-9G(0, 0; x) + 18G(0, 1; x) + 18G(1, 0; x) - 36G(1, 1; x 

-34C(3) + 50C(2)G(-1; y) + 39C(2)G(0; x) + 66C(2)G(1; x) 



96(1 - x)2 

-72C(2)G (-l/y; x) + 24G(0; x)G(-l, -1; y) + 24G(1; x)G(-l, -1; y) 

-36G (-l/y; x) G(-l, -1; y) - 36G(-y; x)G(-l, -1; y) - 18G(0; x)G(0, -1; y) 

+18G (-l/y; x) G(0, -1; y) + 18G(-y; x)G(0, -1; y) + 36G(-1; y)G(l, 0; x) 
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-72G(-1; y)G(l, 1; x) - 18G(-1; y)G (-l/y, 0; x) + 36G(-1; y)G (-1/y, 1; x) 
-18G(-1; y)G(-y, 0; x) + 36G(-1; y)G{-y, 1; x) + 64G(-1, -1, -1; y) 
-36G(-1, 0, -1; y) - 36G(0, -1, -1; y) + 18G(0, 0, -1; y) - 9G(0, 0, 0; x) 
+18G(0, 0, 1; x) + 36G(0, 1, 0; x) - 72G(0, 1, 1; x) + 54G(1, 0, 0; x) 
-108G(1, 0, 1; x) - 108G(1, 1, 0; x) + 216G(1, 1, 1; x) - 18G (-1/y, 0, 0; x) 
+36G (-1/y, 0, 1; x) + 18G (-1/y, 1, 0; x) - 36G (-1/y, 1, 1; x) 



+18G(-2/,l,0;x) -36G(-y,l,l;x) 



(C.22) 



The topology shown in Fig. ^(d) has two Mis. The first one is 




P(i{k'2} Po{ki-k2) Po{ki-pi) Po{ki-pi-p2) Pmih-Ps) Pmih-Ps) 



, (C.23) 



with 




(C.24) 



A-2 



32(1 -x2)(l + y) 

X 



-4C(2) -G(0,0;x) +2G(0,l;x) 

-5C(3) + 8C(2)G(-1; x) + 8C(2)G(-1; y) - 3C(2)G(0; x) 



32(1 -x2)(l + y) 
+16C(2)G(1; x) - 8C(2) G (-1/y; x) - 4G (-1/y; x) G(-l, -1; y) 
+4G(-y; x)G(-l, -1; y) + 2G(0; x)G(0, -1; y) + 2G (-1/y; x) G(0, -1; y) 
-2G{-y- x)G(0, -1; y) + 2G(-1; y)G(0, 0; x) - 4G(-1; y)G(0, 1; x) 
-2G(-1; y) G (-1/y, 0; x) + 4G(-1; y)G (-1/y, 1; x) + 2G(-1; y)G(-y, 0; x) 
-4G(-1; y)G(-y, 1; x) + 2G(-1, 0, 0; x) - 4G(-1, 0, 1; x) + 4G(0, -1, -1; y) 
-2G(0, 0, -1; y) - 3G(0, 0, 0; x) + 6G(0, 0, 1; x) + 2G(0, 1, 0; x) 
-4G(0, 1, 1; x) + 4G(1, 0, 0; x) - 8G(1, 0, 1; x) - 2G (-1/y, 0, 0; x) 
+4G (-1/y, 0, 1; x) + 2G (-1/y, 1, 0; x) - 4G (-1/y, 1, 1; x) - 2G(-y, 1,0; x) 
+4G(-y,l,l;x)] . (C.25) 

The second MI for topology ^(d) is 




^'^ki^'^k2{k2 - Pi - P2? 



Po {^2) Po{ki- k2) Pq {ki -pi) Po{ki -pi -p2) Pm (ki -ps) Pm {k2 -Ps) 



,(C.26) 
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(C.27) 



32(1 + y) ' 
1 

16(1 + y) 
1 



1 



16(l + y)(x + l) L 
+2(x + l)G(-l,-l;y) + (2: 



G(-l;y) 

5C(2)x + 2x- 3C(2) + 2 - 2(x + 1)G(-1; y) 
l)G(0,0;x) -2(x- 1)G(0, 1;2;) 



1 



8(x + 1)(x-1)2g(1;x)G(-1,-1;2/) 



16(1 + y){x + 1)(1 -x + x'^ + xy)l 
+A{x + l)(x - l)2G(-l;y)G(l,0;x) - 8(x + l)(x - 1)2G(-1; y)G(l, l;x) 
-4(x + l)(x - 1)2G(1, 1,0;x) + 8(x + l)(x - 1)^G(1, 1, l;x) 
-8(x2+ya;-x+l) (x - 1)C(2)G(-1; x) + 3 (Sx^+yx -x-l) (x - 1)C(2)G(0; x) 
-4 (x2+4yx-4x+7) (x - 1)C(2)G(1; x) + 8(2/x-x+2)(x - 1)C(2)G (-l/y; x) 
+4(yx - X + 2)(x - 1)G (-l/y; x) G(-l, -1; y) - 4x(2x + y - 1) x 
X (x - l)G(-y; x)G(-l, -1; y) - 2x(2x + y - l)(x - 1)G(0; x)G(0, -1; y) 
-2{yx - X + 2)(x - 1)G (-1/y; x) G(0, -1; y) + 2x(2x + y - l)(x - 1) x 
xG(-y;x)G(0,-l;y) + 2 (x^+yx-x+l) (x - 1)G(0, 0; x) - 2 (x^+yx-x+l) x 
x(x - l)G(-l;y)G(0,0;x) - 4 (x^ + yx - x + l) (x - 1)G(0, 1; x) 
+4 (x^ + yx - X + 1) (x - 1)G(-1; y)G(0, 1; x) + 2(yx - x + 2)(x - 1)G(-1; y) 
G (-1/y, 0; x) - 4(yx - x + 2)(x - 1)G(-1; y)G (-1/y, 1; x) 
-2x(2x + y - l)(x - 1)G(-1; y)G(-y, 0; x) + 4x(2x + y - l)(x - 1)G(-1; y) x 
xG(-y, 1; x) - 2 (x^ + yx - X + l) (x - 1)G(-1, 0, 0; x) + 4 (x^ + yx - x + l) x 
X (x - 1)G(-1, 0, 1; x) + 2x(2x + y - l)(x - 1)G(0, 0, -1; y) + (Sx^ + 3yx 
-3x + l)(x - 1)G(0, 0, 0; x) - 2 (Sx^ + 3yx - 3x + l) (x - 1)G(0, 0, 1; x) 
-2 (x^ + yx - X + 1) (x - 1)G(0, 1, 0; x) + 4 (x^ + yx - X + l) (x - 1)G(0, 1, 1; x) 
-2 (x^ + 2yx - 2x + 3) (x - 1)G(1, 0, 0; x) + 4{x^ + 2yx - 2x 
+3)(x - 1)G(1, 0, 1; x) + 2(yx - x + 2)G (-1/y, 0, 0; x) (x - 1) - 4(yx - x 
+2)(x - 1)G (-1/y, 0, 1; x) - 2(yx - x + 2)(x - 1)G (-1/y, 1, 0; x) + 4(yx - x 
+2)(x - l)G(-l/y,l,l;x) +2x(2x + y - l)(x - l)G(-y, 1, 0; x) -4x(2x + y 
-l)(x - l)G(-y, 1, 1; x) + 2(5C(2)x^ + 4C(3)x3 + 2x^ + 2yx2 + 5yC(2)x2 
-8C(2)x2 + 5yC(3)x2- 4C(3)x2+ 2yx - 3yC(2)x + 8C(2)x + 6C(3)x - 3C(2) - C(3) 



+2) - 2(5C(2)x^ + 2x^ + 2yx^ + AyC,{2)x^ - %{2)x^ + 2yx - AyC{2)x + 7C(2)a; 
-3C(2) + 2)G(-1; y) + 4(x + I) {x'^ + yx - x + l) G(-l, -1; y) - 4(x + l){x^ 
+yx -x + 1)G(-1, -1, -1; y) + 2x{x + l)(y + 1)G(-1, 0, -1; y) - 2(5x3 + 2,yx^ 



-Qx^ -yx + 2x + 1)G(0, -1, -1; y) 



The topology in Fig. ^-(e) involves two Mis. One of them is 



(C.28) 



Poik2)Po{ki -k2)Poiki -pi)Poik2 -Pi)Po{ki -pi -P2)Pm{ki -p^Pm{k2 - 



-,(C.29) 



where 



(C.30) 



j=-4 



A-2 

A. 



12(1 - xY{l + yf 

X 

96(1 - xf{l + yf 

X 



-14G(-1; y) + 9G(0; x) - 18G(1; x) 
48(1 - x)2(l + y)2 L - 12G(-1; y)G(0; x) + 24G(-1; y)G(l; x) 

+4G(-1, -1; y) + 3G(0, 0; x) - 6G(0, 1; x) - 6G(1, 0; x) + 12G(1, 1; x) 

^Sil-xm + yf [-^^^^^) + 46C(2)G(-1; y) - 21C(2)G(0; x) 
+66C(2)G(1; x) - 24C(2)G (-1/y; x) + 24G(0; x)G(-l, -1; y) 
-24G(1; x)G(-l, -1; y) - 12G (-1/y; x) G(-l, -1; y) 
-12G(-y; x)G(-l, -1; y) - 6G(0; x)G(0, -1; y) + 6G (-1/y; x) G(0, -1; y) 
+6G(-y; x)G(0, -1; y) - 12G(-1; y)G(0, 0; x) + 24G(-1; y)G(0, 1; x) 
+36G(-1; y)G(l, 0; x) - 72G(-1; y)G(l, 1; x) - 6G(-1; y)G (-1/y, 0; x) 
+12G(-1; y)G (-1/y, 1; x) - 6G(-1; y)G(-y, 0; x) + 12G(-1; y)G(-y, 1; x) 
+32G(-l,-l,-l;y) -12G(-l,0,-l;y) -12G(0,-l,-l;y) 
+6G(0, 0, -1; y) + 6G(1, 0, 0; x) - 12G(1, 0, 1; x) - 12G(1, 1,0; x) 
+24G(1, 1, 1; x) - 6G (-1/y, 0, 0; x) + 12G (-1/y, 0, 1; x) + 6G (-1/y, 1, 0; x) 
-12G (-1/y, 1, 1; x) + 6G(-y, 1,0; x) - 12G(-y, 1, 1; x)] . (C.31) 

The second MI for the topology in Fig. |3|-(e) is 



Pi- 



D'^A:i2)'^/c2 kl 



-Po (^2) -Po (ki-k-i) Pq (ki-pi) Pq (fc2-Pi) Po {ki-pi-p2) Pm (ki-ps) Pm ik2-P3) 



(C.32) 
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where 



1 -1 



(C.33) 



j=-4 



A. 



A. 



-10G(-l;y) + 3G(0;x) -6G(l;x) , 

-47C(2) - 24G(-1; y)G(0; x) + 48G(-1; y)G(l; x) + 32G(-1, -1; y) 
-85C(3) + 188C(2)G(-1; y) + 96C(2)G(1; x) - 96C(2) G (-1/y; x) 



1 

" 24(1 +y)2 
_ 1 
~ 96(1 + y)2 
1 

~ 192(1 +y)2 

-6G(0,-l;y) 
1 

~ 192(1 +y)2 

+96G(0; x)G(-l, -1; y) - 96G(1; x)G(-l, -1; y) - 48G (-1/y; x) G(-l, -1; y) 
-48G(-y; x)G(-l, -1; y) - 24G(0; x)G(0, -1; y) + 24G (-1/y; x) G(0, -1; y) 
+24G(-y; x)G(0, -1; y) + 48G(-1; y)G(l, 0; x) - 96G(-1; y)G(l, 1; x) 
-24G(-1; y)G (-1/y, 0; x) + 48G(-1; y)G (-1/y, 1; x) - 24G(-1; y)G(-y, 0; x) 
+48G(-1; y)G(-y, 1; x) + 64G(-1, -1, -1; y) - 24G(-1, 0, -1; y) 
-12G(0, -1, -1; y) + 6G(0, 0, -1; y) + 24G(1, 0, 0; x) - 48G(1, 0, 1; x) 
-48G(1, 1, 0; x) + 96G(1, 1, 1; x) - 24G (-1/y, 0, 0; x) + 48G (-1/y, 0, 1; x) 
+24G (-l/y,l, 0; x) -48G (-l/y,l,l; x) + 24G(-y,l, 0; x) -48G(-y, 1, 1; x)] .(C.34) 

FinaUy, the topology in Fig. ^(f) involves three Mis. The first MI we chose is 



Po (fci)Po (A:2)Po (A:i- {k2-pi)Po {ki-pi-p2)Po {k2-pi-p^Prn {ki-pz) 



, (C.35) 



where 



(C.36) 



j=-4 



A. 



24(l-x)4(l+y) ' 
X^ 

96(l-x)4(l+y) 



-10G(-1; y) + 3G(0; x) - 6G(1; x) 



-24- 



A. 



X 



48(l-x)4(l+y) 
,.2 



X 



-5C(2) - 6G(-1; y)G(0; x) + 12G(-l; y)G{l- x)+8G(-l -1; y) 



-13C(3) + 38C(2)G(-1; y) + 9C(2)G(0; x) + 6C(2)G(1; x) 



48(l-a;)4(l+y) 

-24C(2)G (-1/y; x) + 24G(0; x)G(-l, -1; y) - 24G(1; x)G(-l, -1; y) 
-12G(-l/y;x)G(-l,-l;y)-12G(-y;x)G(-l,-l;y) -6G(0;x)G(0,-l;y) 
+6G (-1/y; x) G(0, -1; y) + QG{-y- x)G(0, -1; y) + 12G(-1; y)G(l, 0; x) 
-24G(-1; y)G(l, 1; x) - 6G(-1; y)G (-1/y, 0; x) + 12G(-1; y)G (-1/y, 1; x) 
-6G(-l;y)G(-y,0;x) +12G(-l;y)G(-y,l;x) + 16G(-l,-l,-l;y) 
-12G(-l,0,-l;y) -12G(0,-l,-l;y) +6G(0,0,-l;y)+6G(l,0,0;x) 
-12G(1, 0, 1; x) - 12G(1, 1, 0; x) + 24G(1, 1, 1; x) - 6G (-1/y, 0, 0; x) 
+12G (-1/y, 0, 1; x) + 6G (-1/y, 1, 0; x) - 12G (-1/y, 1, 1; x) + 6G(-y, 1, 0; x) 



-12G(-y,l,l;x) . 
The second MI for the topology Fig. ^(f) was chosen as follows 



(C.37) 



T)'^kiTi'^k2 {ki -piY 



Po (ki) Po ih ) Po ih-k2) Po {k2-Pi) Po (ki-pi-p2) Po (^2-^1-^2) Pm (ki-ps) 



, (C.38) 



where 



Pi = — + (C.39) 

i=-2 



A. 



X 



16(1 - x)3(l + x) 

^2 



4C(2) +G(0,0;x) -2G(0,l;x) 



X 



5C(3) + 24C(2)G(-1; x) - C(2)G(0; x) - 8 C(2)G(1; x) 
- 4G {-l/y; x) G(-l, -1; y) + 4G(-y; x)G(-l, -1; y) 



16(1 - x)3(l + x) 
-8C(2)G(-l/y;x) 

+2G(0; x)G(0, -1; y) + 2G (-l/y; x) G(0, -1; y) - 2G(-y; x)G(0, -1; y) 

-2G(-1; y)G (-l/y, 0; x) + 4G(-1; y)G (-l/y, 1; x) + 2G(-1; y)G(-y, 0; x) 

-4G(-1; y)G(-y, 1; x) + 6G(-1, 0, 0; x) - 12G(-1, 0, 1; x) + 4G(0, -1, -1; y) 

-2G(0, 0, -1; y) + 2G(0, 0, 0; x) - 4G(0, 0, 1; x) - 6G(0, 1, 0; x) + 12G(0, 1, 1; x) 

-2G(1, 0, 0; x) + 4G(1, 0, 1; x) - 2G (-l/y, 0, 0; x) + 4G (-l/y, 0, 1; x) 

+2G (-l/y, 1,0; x) -4G (-l/y, 1, 1; x) -2G(-y, 1,0; x)+4G(-y, 1, 1; x)] . (C.40) 

The last MI for topology |-(f) is 



-25- 



with 



Po {h)Po ik2)Po {ki-k^Po (k2-pi)Po {ki-pi-p2)Po {k2-Pi-P2)Pm (ki-ps) ' 



(C.41) 



m 



4 

i=-A 



(C.42) 



A-2 

A-i 



-8G(-1; y) + 3G(0; x) - 6G(1; x) 



-5C(2) - 3G(-1; y)G(0; x) + 6G(-1; x) + 2G(-1, -1; y) 



-29C(3) + 58C(2)G(-1; y) + 9C(2)G(0; x) + 6C(2) G(l; x) 



192(1 - x)4 

x^ 

96(1 -x)4 
x^ 

24(1 -x)4 L 
x^ 

48(1 - xf 

-24C(2)G (-1/y; x) + 24G(0; x)G(-l, -1; y) - 24G(1; x)G(-l, -1; y) 
-12G {-l/y- x) G(-l, -1; y) - 12G(-y; x)G(-l, -1; y) - 6G(0; x)G(0, -1; y) 
+6G (-l/y; x) G(0, -1; y) + 6G(-y; x)G(0, -1; y) + 12G(-1; y)G(l, 0; x) 
-24G(-l;y)G(l,l;x) -6G(-l;y)G(-l/y,0;x) + 12G(-l;y)G(-l/y,l;x) 
-6G(-l;y)G(-y,0;x) + 12G(-l;y)G(-2/,l;x) +32G(-l,-l,-l;y) 
-12G(-1, 0, -1; y) - 12G(0, -1, -1; y) + 6G(0, 0, -1; y) + 6G(1, 0, 0; x) 
-12G(1, 0, 1; x) - 12G(1, 1, 0; x) + 24G(1, 1, 1; x) - 6G (-l/y, 0, 0; x) 
+12G (-l/y, 0, 1; x) + 6G (-l/y, 1,0; x) - 12G (-l/y, 1, 1; x) + 6G(-y, 1, 0; x) 



-12G(-y,l,l;x) 



(C.43) 
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